Self-organized patterns of macroscopic quantum tunneling in molecular magnets 
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We study low temperature resonant spin tunneling in molecular magnets induced by a field 
sweep with account of dipole-dipole interactions. Numerical simulations uncovered formation of 
self-organized patterns of the magnetization and of the ensuing dipolar field that provide resonant 
condition inside a finite volume of the crystal. This effect is robust with respect to disorder and 
should be relevant to the dynamics of the magnetization steps observed in molecular magnets. 
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Molecular magnets had become focus of interest after 
it was demonstrated [lj that a crystal of Mni2-Acetate 
molecules is equivalent to a regular array of single-domain 
superparamagnetic particles of spin 10. The subsequent 
discovery of the stepwise hysteresis curve in Mni2-Ac @] 
made it a model system for the study of macroscopic 
quantum tunneling of the magnetic moment Q . In a typ- 
ical macroscopic experiment the magnetic field is swept 
at a constant rate, making the spin energy levels |m) 
and |m'} cross at certain fields B = Bk, see Fig. [1] 
Probability of a quantum transition between the cross- 
ing levels affects the size of the magnetization step. In 
this Letter we develop theoretical framework for under- 
standing the dynamics of the steps. We will be concerned 
with the regime of pure quantum tunneling or thermally 
assisted quantum tunneling 0]. 

It had been recognized in the past [f| that the knowl- 
edge of the probability of spin tunneling alone is insuf- 
ficient to describe the tunneling dynamics of molecular 
magnets. This is because the reversal of any individual 
spin changes the long-range dipolar magnetic field of the 
entire crystal, which tunes other spins in or out of res- 
onance. The tunneling dynamics of a molecular magnet 
near a resonance is, therefore, non-local and non-linear 
@. Prokof'ev and Stamp Q suggested that the relax- 
ation could occur near some surface inside the sample 
which is the locus of points characterized by the zero bias 
between spin energy levels. Their local model, however, 
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could not capture the evolution of such a surface in time 
that should occur during magnetic relaxation [7J. Later 
studies of the collective relaxation in molecular magnets 
did take the non-locality into account. They did not 
consider, however, the dynamics generated by the field 
sweep that is important for the formation of a coherent 
spatial structure in the macroscopic tunneling dynamics 
driven by dipolar fields. 

In this Letter we derive and solve non-linear integro- 
differential equations that govern the self-organized tun- 
neling dynamics of spins coupled by dipolar interactions. 

Let n m be the normalized population of the m-th level, 
J2 m Um = 1- I n a typical macroscopic experiment one 
magnetizes the crystal such that initially n_s = 1 while 
populations of other states are zero. The field is then 
swept in the opposite direction at a constant rate. It has 
been demonstrated theoretically [lfj| and confirmed by 
experiments that in the thermally assisted regime 
a step in the magnetization curve is dominated by a 
resonance between particular m and m! levels, see Fig. 
([1]). This resonance is determined by the maximum of 
A 2 mm , exp[-(E m -E- S )/(k B T)], where A mm , is the tun- 
nel splitting of the m-th and m'-th levels, E m is the en- 
ergy of the m-th level, and T is temperature. On ap- 
proaching this resonance, the population of the | — S) 
state satisfies the following equation [1] 
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Here W is the energy bias between m and m', and 
T m , m > = T m + r m /, T m being the rate of transitions from 
to to the neighboring spin levels. The overdamped condi- 
tion T m . m i 3> A is assumed, which is the case of practical 
interest. Since the m'} state quickly decays down to the 
\S) state, the magnetization is given by 



FIG. 1: Thermally assisted tunneling. 



<J Z EE (S,)/S = 1 - 277_5 . (3) 
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At a constant sweep rate, W = Vwt, when T m ,m' is small 
compared to the distance between adjacent spin levels, 
one can integrate Eq. |T]) from — oo to oo to obtain the 
population of the | — S) state after the level crossing. 
This gives the Landau-Zener like asymptotic result, 

n_ s (oo) = e~ e , (4) 

with e corrected for thermally assisted tunneling 
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The first magnetization step occurs when the tunnel 
splitting of the two states is sufficient to provide a not 
very small e for a given sweep rate vyy. All of the above 
formulas apply to the case of pure quantum tunneling 
that occurs at T = 0. In the latter case m — ~S. 

The above simple picture contradicts experiments. In- 
deed, according to Eq. ([2]) the width of the magnetiza- 
tion step due to the relaxation dynamics described by 
Eq. ([1]) is determined by T m > . Such a width is very small 
compared to the experimental value. One could naively 
think that the relaxation of any individual molecule is 
still determined by Eqs. (JT|) and @ but the resonant 
field is spread due to inhomogeneous dipolar and hyper- 
fine fields, and due to structural disorder. This, however, 
cannot explain the height of the magnetization step as 
this would typically result in the total magnetization re- 
versal during the first step. As we shall see, the non- 
trivial space-time dependence of the local energy bias W 
due the evolution of the local dipolar field B( D ) solves 
this controversy. 

To simplify the problem we stick to the geometry of a 
long cylinder of length L and radius R <C L, and restrict 
our consideration by the longitudinal components of the 
fields. The transverse component of the field enters the 
problem through its known effect on A mm / . The bias at 
the cite i is given by 



W t = (to' - to) 5Mb ( B - B k + Bl B) ) ee W cxt + W t 
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where Bk is the resonant value of the external field for 
the fc-th tunneling resonance (fc = S — m') in the absence 
of the dipolar field. The dipolar component of the bias 
is given by 



W * (D) = ( 1 " 7n~ ) Ed A ' A = £ hi* 3 ,, (7) 



where Ey> = (g/iB'S) /vq is the dipolar energy, iib/fcB 
0.0671 K for Mni2-Ac, vq is the unit-cell volume, and 
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is the dimensionless dipole-dipole interaction between the 
spins at the sites i and j ^ i. 
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FIG. 2: Color online: Moving wall of resonant spin tunneling 
induced by a slow sweep, e = 10. There are well-developed 
spatially quasiperiodic, time independent structures behind 
the wall. 



To write down the dynamical equation for the entire 
sample, it is convenient to use dimensionless variables 
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is the relaxation rate on resonance, see Eq. j2| . Equation 
(fTj) becomes a non-linear integro-differential equation, 



dt 



n-s(z, t) = —F(z, t)n~s(z, t) . 



(11) 



Here F contains integral dependence on n_g{z 1 1) via W, 
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Summation of contributions from all sites in Eq. J7J gives 
the following expression for a continuous variable W(z, t): 



W = W cxt + W {D ^ = W cxt - k[1 - 2n_ s (z, t)] 



L/(2R) 



(? -Z) +1 



3/2 



(13) 



where v is the number of molecules per unit cell and 

« = -5--2^^i- ( 14 ) 



The summation in the last formula goes over the volume 
of a large sphere. For the body-centered tetragonal lat- 
tice of Mni2-Ac one has v — 2 and fc = 14.6 fl3| . In 
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FIG. 3: Color online: Spatial profiles of the metastable pop- 
ulation n-s and the reduced bias W at Wcxt = 4.3 and 
-Ed = 200. Everywhere in the wall the system is near the 
resonance, W ~ 0. 



the depth of a uniformly polarized elongated sample Eq. 
(flU)) provides the dipolar field of 526G in good agree- 
ment with the measured value of 515 ± 85G 
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tions (|li p -(|13 p should be solved with the initial condition 
ns(t = — oo) = 1. 

For Mni 2 -Ac Ejy is a large parameter, En > 100, so 
that one can naively think that F is always small except 
for a very brief period of time when the field sweep brings 
the bulk of the sample close to the resonance. If this was 
true, the total relaxation would have been vanishingly 
small. However, as we shall see below, the system finds 
the way to relax faster by forming a moving wall of fi- 
nite width I ~ R inside which W is so small that F is 
not significantly reduced from its maximal value F = 1. 
Inside the wall region the resonant tunneling transitions 
take place. The greater is Ejy, the closer to zero is the 
energy bias W in the wall. Beyond the wall region W 
deviates strongly from zero, F becomes very small, and 
the relaxation effectively freezes. 

Numerical solution of Eq. ([Til based upon its dis- 
cretization over the length of the sample yields a propa- 
gating wall of tunneling shown in Fig. [2j Spatial depen- 
dence of the metastable population n_s and of the energy 
bias are shown in Fig. [3l Striking universal features of 
the relaxation process have been uncovered by the simu- 
lations. The wall ignites as the bias reaches the "magic" 
value of W C xt = 4.3, which corresponds to B = £?fc + 194G 
in Mni2-Ac regardless of the sweep rate and the level m 
that dominates the transition. One can see from Fig. [3] 
that everywhere in the wall the system is near the reso- 
nance, W w 0. The reduced speed of the wall, 

(15) 



-nT rcs 

depends on the parameter Ed . The propagation of the 
wall does not result in the total magnetization reversal. 
The metastable population behind the wall initiated by 
a slow sweep is a universal number rif « 0.32, see Fig. 3] 




Relaxation of the average 
metastable population 
in the sample under a slow 
sweep of the external bias 



FIG. 4: Time dependence of the average metastable popula- 
tion fi-s = (n-s)- 




z/(2R) 



FIG. 5: Color online: Spatial profiles of the metastable pop- 
ulation n-s and the reduced bias W in the presence of a 
relatively strong disorder. The condition W ~ inside the 
wall is not affected by disorder. 



The width of the wall, 21, is of the order of the diameter 
of the elongated sample. Figs. [5] and show the tunneling 
wall and the quasiperiodic structure behind the wall in 
the left half of the sample. In a totally symmetric sample, 
two walls will simultaneously ignite on the left and on 
the right sides of the sample, and then move towards 
the center. Deviation from the dynamics depicted in the 
above figures will occur only when the distance between 
the two walls becomes comparable to R. For L ^> R this 
should not change our conclusions. 

So far we have not considered the effect of disorder 
on the formation of the tunneling wall. For Mni2-Ac 
this question is important because nuclear spins, solvent 
disorder, and crystal defects produce randomness in the 
local values of Bk and local dipolar fields. To study the 
effect of disorder we added a random component, Wrand) 
to the energy bias at each lattice site. While disorder 
produces visible fluctuations of n_s and W behind the 
wall and ahead of the wall, the existence of the wall with 
W « is not affected by disorder, see Fig. 

We believe that our results are relevant to the magne- 
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tization steps commonly observed in molecular magnets 
at a slow field sweep. These results are robust with re- 
spect to other activation scenarios, e.g. those involving 
manifolds with a different total spin S [l5| . The latter 
will simply modify the expression for r rcs , Eq. p0[) . 

Without dipolar fields the energy width of the magne- 
tization step, according to Eq. ((2]), is AWr = hT m . m i/2. 
With account of the dipolar fields it is determined by 
the time it takes the tunneling wall to cross the sample, 
t = L/v. During that time the energy bias changes by 
AW — vy/t = vwL/v = vw L / (v* RT TCS ) . This gives the 
following ratio of the field widths of the magnetization 
step with and without the dipolar fields: 



AB 

ABr 



TT 

V*£ 



where we have used Eqs. ([5]) and (fT0|) . For the parameters 
of Mni2-Ac, in the limit of a very slow sweep, one numer- 
ically obtains v* ~ 1CT 3 . At e = 30 and L/R = 30 used 
in the simulations, Eq. (fTtj)) then gives AB/ABr ~ 10 3 . 
This brings the width of the step in the range that is 
usually observed in experiment. So far we have studied 
the dynamics of the first magnetization step. The second 
step begins with a quasiperiodic structure of magnetiza- 
tion left after the first step. This dynamics will be studied 
elsewhere based upon the equations derived above. 

After this work was completed, it had come to our at- 
tention that the oscillating magnetization near tunneling 
resonance, depicted in Fig. [3l has been already seen in 
experiment [161 ]. We encourage experimentalists to per- 
form detailed local measurements of molecular magnets 
to confirm that the magnetization step involves propa- 
gating walls of quantum spin tunneling. 

This work has been supported by the NSF Grant No. 
DMR-0703639. 
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